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Hong-Ou-Mandel (HOM) interference between independent photon sources (HOMI-IPS) is the 
fundamental block for quantum information processing, such as quantum gate, Shor’s algorithm. 

Boson sampling, etc. All the previous HOMI-IPS experiments were carried out in time-domain, 
however, the spectral information during the interference was lost, due to technical difficulties. Here, 
we investigate the HOMI-IPS in spectral domain using the recently developed fast fiber spectrometer, 
and demonstrate the spectral distribution during the HOM interference between two heralded single¬ 
photon sources, and two thermal sources. This experiment can not only deepen our understanding 
of HOMI-IPS in the spectral domain, but also be utilized to improve the visibility by post-processing 
spectral filtering. 


I. INTRODUCTION 

Hong-Ou-Mandel (HOM) interference [1] between inde¬ 
pendent photon sources (HOMI-IPS) is at the heart of 
quantum information processing involving the quantum 
interference of single photons laig. Scaling such sys¬ 
tems to large dimensions requires many single-photons 
to interfere in a single mode H E]. This establishes a 
fundamental challenge since the outcome fidelity of the 
quantum processor depends on the interference visibility 
of all photons in all modes. To insure high output-fidelity 
the interference of at least two independent photons must 
approach unity. 

A variety of HOM-IPS experiments have been demon¬ 
strated at near-infrared wavelengths [SHIS] or telecom 
wavelengths pMll2Q] . and with heralded single-photon 
sources jSa IIBHIDI, weak coherent sources mm , or 
even thermal sources ns mi- In all the previous HOMI- 
IPS experiments, the HOM dips are obtained by record¬ 
ing the coincidence counts as a function of the temporal 
delay. The knowledge of the spectral correlations of the 
interfering photons can reveal some important informa¬ 
tion about the interference visibility [22]. In the past, the 
spectral correlation content could not be observed due to 
the long acquisition times required when measuring spec¬ 
tral correlations using two tunable bandpass filters. 

Recently, a new technique for measuring spectral cor¬ 
relations in a HOM interference was presented by Ger- 
rits, et al |2S||2l]. This technique has the merits of high 
speed and high spectral resolution, and therefore allows 
measuring and analyzing the spectral correlations from a 
HOM interference. 

Here, we measure the spectral correlation between 
the two output ports of the HOMI-IPS beam splitter 
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with two different photon sources (heralded single-photon 
sources, and thermal sources) using the method devel¬ 
oped in [23]. This experiment can help deepen our un¬ 
derstanding of HOMI-IPS from the viewpoint of spectral 
domain and allows improving the HOM-IPS visibility by 
post-processing spectral filtering. 

First, we present the theory and simulation for HOM 
interference between two heralded single-photon states, 
and between two thermal states. We then present our 
experimental results for these two kinds of HOMI-IPS ex¬ 
periments. The appendix of this manuscript contains our 
theoretical model for HOM interference between two her¬ 
alded single-photon states, by considering the infinitely 
higher order components in the photon sources. 


II. THEORY AND SIMULATION 

A. Theory of four-fold HOM interference between 
two heralded single-photon states 

Two independent spontaneous parametric down con¬ 
version (SPDC) sources are used in the experiment. The 
signal and idler photons from the first and the sec¬ 
ond SPDG source have the joint spectral amplitude of 
fi{uJsi^^ii) and f 2 {^s 2 ^^i 2 )^ where uj is the angular fre¬ 
quency. The subscripts 1 and 2 represent the first and 
the second SPDG, respectively; s and i denote the down- 
converted photons, i.e., the signal and idler. The signal 
photons {si and 52) are sent to a 50/50 beamsplitter for 
HOM interference, while the idler photons (il and i2) 
function as heralders, i.e., detected by single-photon de¬ 
tectors to claim the existence of 5 I and s2 (see the layout 
in Fig. 1^. The four-fold coincidence probability P 4 (r) 
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as a function of the time delay r can be written as [22] : 

oo oo oo oo 

P4(T) = i //// du}sidu}s2duJiidiOi2h{T), (1) 


0 0 0 0 


where, 


hir) = |/i(wsi,a;ii)/2(ws2,Wi2)- 

/l(Ws2,Wii)/2(Wsl,Wi2)e“*(“"2-a;si)T|2_ 


( 2 ) 


The correlated spectral intensity (CSI), / 4 (t), for the two 
beam splitter output ports at time delay r can be ex¬ 
pressed as: 


/ 4 ( t ) = J J duJiiduJi2h{r). 

0 0 


(3) 


Figure shows the simulated HOM dip and CSI us¬ 
ing the experimental parameters of the setup: ~ 2 ps 
pump pulses at 792 nm, 30-mm-long periodically poled 
KTi 0 P 04 (PPKTP) crystal. As shown in Fig. [^a), 
without any bandpass filtering, the visibility of the four¬ 
fold HOM dip is 82.0%, and the FWHM (full width at 
half maximum) is 4.91 ps (1.47 mm). The intrinsic spec¬ 
tral purity of this source is 0.82 [25l |26], which is also 
the upper bound of the four-fold HOM interference vis¬ 
ibility in our case laiiHi- The theoretical prediction in 
Fig. gb -g) shows fringes in the CSI at different delay 
positions. The fringe frequency increases with increasing 
distance from the center of the HOM dip. 


B. Theory of two-fold HOM interference between 
two thermal states 


where (0) = ^(oo) = 0, and S’ < ^. For simplic¬ 

ity, here we assume two nonlinear crystals are identical: 

= / 2 (Ws 2 ,Wi 2 ) = f{uJs,UJi). is the 

complex conjugate of /(cjs,cj^). In Eq.®, ^ — S{t) 
is the interference parts, which is contributed by single 
photons from two different crystals, and ^ S S is the 
background parts, which is contributed by two photons 
from one crystal. 

In fact, hir) in Eq.([^ can be expressed as m 


= sz/— S{t). (9) 

Eq. (§ and Eq. suggest the HOM interference between 
two thermal states is including the HOM interference be¬ 
tween two heralded single-photon states ^ — S{r), plus 
a constant background term of ^ S S. The correlated 
spectral intensity (CSI), for the two beam splitter 

output ports at time delay r can be expressed as: 


oo oo 

t('r) = y y dcOidoJiJtir). (10) 

0 0 

With the same condition for Eig. we simulate the 
two-fold HOM dip in Eig. [^a), and /t(r) at different 
delay positions are shown in Eig. ib -g). The visibility 
V of HOM interference between two thermal state has 
an upper bound of ^ EH: 


Pi(oo) - PtjO) ^ ^ ^ 1 

Pf(^oc^ S 3 


(11) 


where Pt{±oc) (x + S) and Pt(0) (x ^(2^). Re¬ 

fer to m for more details of the theory of interference 
between two thermal states. 


The two-fold coincidence probability Ptir) between 
signal 1 and signal2 (without the gating of idler 1 and 
idler2, signall and signal2 are thermal states) as a func¬ 
tion of the time delay r can be written as [27] 


III. EXPERIMENT AND RESULTS 

A. Experiment of four-fold HOM interference 
between two heralded single-photon states 


oo oo CX3 oo 


P((r) ~ \ j j j j d^sduJs.diOiduJi.ItiT), 


0 0 0 0 


(4) 


where. 


Ii{t^ — -j- S -j- — S{t^^ (5) 

with 

= 2 |/(ws,Wi)/(ws,,Wj,)y (6) 


S{r) = 

( 8 ) 


The experimental setup is shown in Eig. The ex¬ 
perimental setup is the extension of our previous work 
m- We add two kinds of dispersion fibers to the four 
detection channels. Each output port of the fiber beam 
splitter (EBS) is connected to a dispersion compensa¬ 
tion fiber module (DCEM), a 7.53-km-long single-mode 
fiber (SME) with high dispersion designed to compensate 
the dispersion of a 50-km-long commercial SME at the 
L-band of telecom wavelengths. The DCEM has a dis¬ 
persion of 125.0 ps/km/nm and an insertion loss of 5.2 
dB. Due to the large dispersion, the DCEM converts the 
spectral distribution information of the photons into ar¬ 
rival time information, which can be determined with 
high accuracy using superconducting nanowire single¬ 
photon detectors (SNSPDs). The DCEMs, SNSPDs 
and the time interval analyzer (TIA) constitute a fiber 
spectrometer [3T1 [32] . To ensure that the idler photons 
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FIG. 1. The simulation for the HOM interference between two heralded single-photon states.(a) The simulated HOM dip. 
(b-e) Simulated correlated spectral intensity (CSI) at different delay positions in the HOM dip. 
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FIG. 2. The simulation for the HOM interference between two thermal states, (a) The simulated two-fold HOM dip.(b-g) GSI 
at different delay positions in the HOM dip. 


have the same optical path delay as the signal photons, 
we add two 7.53-km-long commercial SMFs (L-SMF in 
Fig.|3| to the idler channels. L-SMF has a dispersion of 
27.3 ps/km/nm and an insertion loss of 1.6 dB. Assum¬ 
ing a system jitter of 100 ps, the resolution of our fiber 
spectrometer is approximately 0.11 nm. 

We pump the PPKTP crystal sources with a 
TiiSapphire oscillator at 76 MHz, each with an average 
pump power of 100 mW. We obtain a singles count rate 
of ^ 1.6 Mcps, a two-fold coincidence count rate (2f-CC) 
of 350 kcps, and a four-fold coincidence count rate (4f- 
CC) of 880 cps. To our knowledge, this 4f-CC is the 
highest rate reported at telecom wavelengths. The mean 
photon number per pulse and the total efficiency of each 
spatial mode are then estimated to be n = 0.114 and 
T] = 0.19, respectively. After insertion of the DCFM and 
the L-SMF, the 2f-CC and 4f-CC are reduced to 35 kcps 
and 37.7 cps. First, we measure a HOM interference 
between signall and signal2, with idler 1 and idler2 as 
heralders. The result is shown in Fig. I^a). Without 
subtraction of any background counts, the raw visibil¬ 
ity of the HOM dip in Fig. [^a) is 68.9±1.3%, which 
is consistent with our previous results in [181 ES] . The 
degradation from the expected visibility of 82.0% (Fig. 


[^a)) to 68.9±1.3% is mainly caused by the multi-pair 
emission in the SPDC [34]. This is confirmed by the 
numerical simulation. Based on the theory in [34], we 
develop a theoretical model including losses, multi-pair 
photon emissions (up to infinitely higher order photons), 
and the realistic frequency mode structure, i.e. the joint 
spectral amplitude (see Appendix). For n = 0.114 and 
T] = 0.19, the visibility is simulated to be 68.1% which 
agrees with the experimental result. Note that in princi¬ 
ple, the visibility can be further improved by decreasing 
the pump power, but with the tradeoff of lower count 
rate. The FWHM of this HOM dip is 4.44 ps (1.33 mm), 
which corresponds well with the theoretical prediction in 
Fig. [^a). 

Figure [^b-g) shows the measurements of the CSI at 
different delay positions, 0 FWHM to 5 FWHMs away 
from the dip center. The fringe observed at the HOM 
dip is very weak (shown in Fig. [^b)), as expected from 
Eq. <§■ The CSI clearly shows fringes at delay positions 
from IFWHM to 4FWHM in Fig. gc-f). The shape of 
the CSI is almost circular at 5FWHM (Fig. Qg)), and 
fringes are not observed due to the limited resolution of 
our fiber spectrometer. The observed CSI from OFWHM 
to 4FWHM agree well with our theoretical predictions in 
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FIG. 3. The experimental setup. Picosecond laser pulses (76 MHz, 792 nm, temporal duration ~ 2 ps) from a mode- 
locked Titanium sapphire laser (Mira900, Coherent Inc.) were divided by a polarizing beam splitter (PBS) into two paths, 
and pumped two 30-mm-long PPKTP crystals with a poling period of 46.1 /xm for type-II SPDC. The downconverted signal 
(horizontal polarization) and idler (vertical polarization) photons with a degenerate wavelength of 1584 nm were divided by 
two PBSs, and then coupled into single-mode fibers (SMFC). The signal photons were sent to a 50/50 fiber beam splitter 
(FBS) and connected to two dispersion compensation fiber modules (DCFM), in order to generate large dispersion. The idler 
photons were connected to two long-distance single-mode fiber spools (L-SMF), to achieve the same path delay as the signal 
photons. Finally, all collected photons were sent to four superconducting nanowire single-photon detectors (SNSPDs), which 
were connected to four input channels of a time interval analyzer (TIA).Our SNSPDs have a system detection efficiency of 
around 70% with a dark count rate less than 1 kcps [281130] . Electrical signals recording the timing information of the pump 
pulses were sent to the TIA for synchronization (Sync). LPF = long pass filters, HWP = half-wave plate, QWP = quarter-wave 
plate, FPC = fiber polarization controller. 
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FIG. 4. (a) The measured four-fold HOM dip between signall and signal2, with idlerl and idler2 as herladers. (b-e) Measured 

correlated spectral intensity (GSI) at different delay positions in the HOM dip. Each GSI data was accumulated in 2 hours. 


Fig. Qb-f). 

B. Experiment of two-fold HOM interference 
between two thermal states 

In the following we present our results when measur¬ 
ing the two-fold coincidences between signall and signal2, 
without the heralding of idlerl and idler2, i.e., two ther¬ 
mal states. The measured HOM dip is shown in Fig. 
I^a), with a visibility of 24.6±0.3% and FWHM of 4.62 
ps (1.39 mm). This obtained visibility is comparable to 
our previous results in m- The measured CSI data are 


shown in Fig. [^b-g), with accumulation time of 10 min¬ 
utes for each data set. Surprisingly, the CSIs also show 
fringes at the different delay positions from IFWHM to 
4FWHM, even with an interference visibility is as low as 
24.6±0.3%. These measured results correspond well with 
our simulated results in Fig. The degradation from 
the expected visibility of 29.1% (Fig. [^a)) to 24.6±0.3% 
is mainly caused by the limited indistinguishability of the 
two different SPDC sources. 
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FIG. 5. (a) The measured two-fold HOM dip between signall and signal2. (b-e) Measured CSI at different delay positions in 

the HOM dip. Each data was accumulated in 10 minutes. 


IV. DISCUSSION AND OUTLOOK 

It is meaningful to compare the HOMI-IPS (using two 
SPDCs, e.g. the case in this experimeut) to the HOM 
interfereuce using a twiu photon source (usiug ouly oue 
SPDC, e.g. the case iu ID)- Iu the latter case, the two¬ 
fold coiucideuce probability p 2 (^) as a fuuctiou of the 
time delay r iu the HOM iuterfereuce is giveu by j35F 

El]: 


c» oo 

P2{t) = \ j J dojsduji [f{u}s,u}i) - 

0 0 

( 12 ) 

The CSI of the two beam splitter output ports /(r) iu 
this HOM iuterfereuce cau be expressed as: 


hir) = \[f{us,Ui) - f{LOi,LOs)e (13) 

To achieve perfect HOM iuterfereuce, Eq.([T^ re¬ 
quires symmetry of the joiut spectral distributiou, i.e., 
f{ujs,uji) = cjs), but doesu’t require high spectral 
purity. However, to achieve high visibility, Eq. 0 re¬ 
quires uot ouly high symmetry, but also a high spec¬ 
tral purity. Iu Eq. 0, to obtaiu 100% visibility, i.e., 
P 4 ( 0 ) = 0, the couditiou of ^ii)/ 2 (^s 2 ,^ 22 ) = 

/i(<^s 2 , ^ii)/ 2 (<^si 5 ^ 22 ) must be satisfied. To achieve this 
couditiou for all cj, firstly the two SPDC sources must 
be equal (iudistiuguishable): fi{00s^00i) = = 

aud secoudly the two sources must be fac¬ 
torable (spectrally pure): = gs{^s)9i{^i) ED- 

Besides the observatiou of the CSI, the techuique iu 
this work cau also be applied to improve the visibility 
of HOMI-IPS by post-processiug spectral filteriug. Here, 
we cau improve the visibility from 68.9±1.3% (iu Fig. 
|4];a)) to 75.7±1.3% (in Fig. ||;a)) and 77.5±1.5% (in 
Fig. §b)), by decreasiug the coiucideuce wiudow from 
5 us to 2 us aud 1 us, correspoudiug to spectral widths 
from 2.66 um to 1.06 um aud 0.53 um, respectively. The 
HOM iuterfereuce visibilities with differeut coiucideuce 
wiudow widths are showu iu Fig. 0c). Fig. 0 c) clearly 


shows that the visibility iucreases as the width of the 
coiucideuce wiudow decreases. Thus, uarrowiug the co¬ 
iucideuce wiudow acts as a spectral filter. Wheu uarrow¬ 
iug the coiucideuce wiudow, the spectral purities of the 
sources are improved aud the HOM visibility is improved. 
Thus, oue cau utilize temporal filteriug to euhauce the 
visibility of the HOM-IPS. However, at the cost of lower 
photou flux. 


V. CONCLUSION 

Iu summary, we have theoretically simulated aud ex- 
perimeutally demoustrated the CSI from a four-fold 
HOM iuterfereuce betweeu two heralded siugle-photou 
sources. We also demoustrated the CSI iu a two-fold 
HOM iuterfereuce betweeu two thermal sources. The 
iuterfereuce patteru iu the CSI is clearly observed aud 
matches our theoretical predictious well. This experi¬ 
meut cau help euable the uuderstaudiug of the HOMI- 
IPS iu the spectral domaiu, aud the techuique cau be 
utilized to improve the visibility of HOMI-IPS by tem¬ 
poral filteriug. 
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APPENDIX: THE THEORETICAL MODEL OF 
THE FOUR-FOLD HOM INTERFERENCE 

This appeudix briefly describe our theoretical model 
of the four-fold HOM iuterfereuce. The model iu- 
cludes losses iu the setup, the multi-pair emissiou (up 
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FIG. 6. (a) Post-processed HOM dip with coincidence window of 2 ns (1 nm). (b) Post-processed HOM dip with coincidence 

window of 1 ns (0.5 nm). (c) Post-processed HOM dip visibility as a function of coincidence window. 


to infinitely higher order), and the frequency multimode 
structure of the pulsed photon-pair source. Note that 
in the theoretical model of heralded single-photon state 
for Eq. 0 and Eq.([T^, only the one-pair component in 
SPDC is considered. While in the theoretical model of 
thermal state for Eq. 0 , both the one-pair and two-pair 
components in SPDC are included. 

An SPDC theory including the first two (losses and 
multi-pair emission) was developed in [3lj where the 
quantum state emitted from the SPDC source is de¬ 
scribed as a Gaussian continuous variable state, more 
precisely two-mode squeezed vacuum, which is dis¬ 
tributed in a single frequency mode. Eor pulsed SPDC 
sources, however, photon-pairs are emitted with more 
complicated frequency distribution and thus has a multi- 
mode structure in frequency. This is particularly im¬ 
portant to simulate the experiments including the inter¬ 
ference between the photons from independent sources 
lau. Here we show how the frequency multimode struc¬ 
ture can be incorporated into the theory in [34] (complete 
description of the theory will appear elsewhere). 

The quantum state emitted from a pulsed SPDC 
source is described by [35] . 


|T) = exp 


•J dus j (kjif{u)s,u)i)al{u)s)d\{u)i) - h.c. 


| 0 ), 


. (14) 

where r is the squeezing parameter (pump power), d|(ci;s) 
and a\{uji) are creation operators for the signal and idler 
modes with frequencies uJq and respectively. 
is the joint spectral amplitude, which is a product of the 
pump distribution and the phase-matching function of 
the nonlinear crystal m ESI Eg. We can assume r to 
be real and positive without losing generality. This joint 
spectral amplitude can be decomposed by the Schmidt 
decomposition as IBEHl 


f{0Js,0Ji) = ^ \/^igi{oJs)hi{uJi), (15) 

i 


where \i is the Schmidt eigenvalue for the Ith Schmidt 
mode. This decomposition allows us to describe the state 


in Eq. (14) as 


where 


= exp r — h.c. 

i 

|0) 


= H®^p ~ 

1 

|0) 


= |«'(rv^))|«'(rv^))--- 
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(16) 

b] = J du}sgi{ujs)dl{u}s), 


(17) 

II 


(18) 


Note that bi and q satisfy a standard bosonic commuta¬ 
tion relation = Sim- In the second line 

of Eq.( 16), decomposition of the exponential term is pos¬ 
sible since the Schmidt modes are orthonormal. The last 
line of Eq.([l^ represents that in the Schmidt mode basis, 
the state is described by a tensor product of two-mode 
squeezed vacuum |T(rV^))‘ 


l^(^vT)) = ^ (tanhryA)') |n)Bjn)c,, 

V L ^ 

(19) 

with the effective squeezing parameter r^/Xl. Therefore, 
to extend the single frequency mode based theory in [34] , 
what need is simply consider a tensor product of multiple 
two-mode squeezed vacua (see Eig. [^. 

Distribution of the Schmidt eigenvalues are theo¬ 
retically calculated using the experimental parameters 
of our setup: ^ 2 ps pump pulses at 792 nm, 30- 
mm-long periodically poled KTiOP 04 (PPKTP) crys¬ 
tal. Under this condition, it can be calculated that 
{Ai,A 2,...} = {0.9,0.025,0.025,0.01,0.009,0.004...} 
[25l|39lllQ]. In our numerical simulation, we include the 
first five Schmidt modes and truncate the rest. 

To fit the simulation with the experiment, we also need 
to connect r with the experimentally observable param¬ 
eter, for example, the photon-pair generation rate of the 






























7 



signal 


idler 


[34] . Instead of a single TMSV in [34] , we start from the 
covariance matrix of two TMSVs and then after applying 
beam splitting operation with mode matching factor 
we obtain the four-spatial-mode covariance matrix 7 (r, 
where r is the squeezing parameter of the TMSVs (see 
[34] for the details of the covariance matrix approach). 
The HOM visibility is defined by 


FIG. 7. (a) single frequency mode SPDC source, (b) pulsed 

(multiple frequency mode) SPDC source in Schmidt mode 
expansion. 


V = 


( 22 ) 


SPDC source. The latter corresponds to the probability 
that the source emit non-zero photons: 


P = 


= l-n|(00|vl/(rv7)) 

I 


( 20 ) 


Pinging Eq. (19) into, Eq. (20), we obtain the relation 

( 21 ) 


p = 1 — ^coshry/X^ 

i 


which allows us to derive r numerically from experimen¬ 
tally estimated p. Note that the photon pair generation 
rate p is related to the mean photon number per pulse n 
via p = h/(h + 1). 

The calculation of the four-fold HOM dip is then car¬ 
ried out by extending the two-fold HOM dip analysis in 


where Pmean and Pmin are the four-fold coincidence count 
probabilities with and without optical path delay, respec¬ 
tively. These coincidence counts are then given by 


P = 1 

mean 


E 


4C*m 


(- 2 )*^ 


m=l C(4,m) ^riidet (7C(4,m)(?’i,0) + J) 

(23) 

Pmin is also obtained by a similar expression with 
7c(4,m)(n,l) instead of 7c(4,m)(?’i,0). 7c(4,m)(?’,e is the 
submatrix of 7 (r, ^) taking modes C(4, m). C(4,m) is 
the m-detector combination from the four detectors (see 
Sec. 3.3 of [34|). = r^/Xi is the effective squeezing pa¬ 

rameter for the ith Schmidt mode where we include the 
first five modes. The HOM visibility is then obtained by 
plugging the above Pmean and Pmin into Eq. (22). 
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